Spin- 1/2 Heisenberg antiferromanget on kagome: a Z 2 spin liquid with fermionic spinons 
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Motivated by recent numerical and experimental studies of the spin- 1/2 Heisenberg antiferromagnet on 
kagome, we formulate a many-body model for fermionic spinons introduced by us earlier [Phys. Rev. Lett. 
103, 187203 (2009)]. The spinons interact with an emergent U(l) gauge field and experience strong short- 
range attraction in the 5 = channel. The ground state of the model is genetically a Zn liquid. We calculate 
the edge of the two-spinon continuum and compare the theory to the slave-fermion approach to the Heisenberg 
model. 



I. INTRODUCTION 

The spin- 1/2 Heisenberg antiferromagnet on kagome 
(Fig. 1) has been extensively studied for over two decades. 12 
Combining geometrical frustration and strong quantum fluctu- 
ation, the model is expected to host unconventional magnetic 
orders as well as fractionalized excitations. Several Cu 2+ 
based kagome materials have been synthesized. Among them, 
one of the most promising realization 3 of the model is the 
herbertsmithite, ZnCuaOHgC^. Bulk susceptibility 3 as well 
as nuclear magnetic resonance (NMR) studies 4 5 of the her- 
bertsmithite conclude that the ground state is a gapless spin 
liquid. Inelastic neutron scattering study 6 has revealed broad 
diffusive structure factor instead of sharp features associated 
with magnons. Relevant perturbations are anisotropic interac- 
tions, especially the Dzyaloshinski-Moryia (DM) term, 7 and 
the presence of paramagnetic impurities. 8 

On the theory side, pioneering two-dimensional density 
matrix renormalization group (DMRG) calculations by sev- 
eral groups 9-12 have established that the ground state of the 
model is a Z 2 spin liquid with a finite spin gap of approxi- 
mately 0.1J, where J is the strength of the exchange coupling. 
Although the finite spin gap seems to contradict experiments, 
it is hoped that the inconsistency can be resolved by taken 
into account the additional perturbations in materials. While 
the effect of DM interaction can certainly be tackled by pow- 
erful numerical techniques such as two dimensional DMRG, 
the effect of Cu substitutions can be tricky to obtain due to the 
size constraint of the finite cluster. 

An analytical framework is thus needed to bridge the gap. 
Ideally, such a theory should yield a Z2 liquid ground state 
and incorporate the effects of lattice defects and other pertur- 
bations. One possible route is offered by slave-particle ap- 
proaches, in which spin variables S = (S x , S v , S z ) are ex- 
pressed in terms of new fictitious particles: 



Q 1 t 

S = 2 a a <Ta P a P> 



(1) 



where the flavor index a — ±1/2 labels particles with up 
and down spins, and <r a p is the triplet of Pauli matrices. For 
either bosonic or fermionic ladder operators a? a and a a , the 
variables (1) satisfy the spin commutation relations. The spin 
length S is set to 1/2 by fixing the net number of particles 



FIG. 1: (a) The kagome lattice with dimers (black thick bonds), a 
defect triangle (the red triangle) and a spinon (the blue spin). The 
corresponding arrow representation is displayed in part (b). 



on a site, a^ a a a = 1. The binary exchange interaction of the 
Heisenberg model translates into quartic interactions between 
the new particles, so that the new problem is not easier to solve 
than the original one. Solutions are usually based on a mean- 
field approximation, in which the constraint on the particle 
numbers is satisfied only on average, (a^a Q ) = 1. The mean- 
field solution is justified if the number of particle flavors N is 
increased from two to infinity. For a finite N, and particularly 
for the physical case N = 2, enforcing the particle-number 
constraint requires going beyond a mean-field treatment. 

For the Heisenberg antiferromagnet on kagome, the 
bosonic route yields a ^-liquid ground state with magnetic 
excitations carrying spin 1/2, known as spinons. Like the 
particles from which spins are constructed, the spin excita- 
tions exhibit the Bose statistics. 13 Adding DM interactions of 
sufficient strength induces a quantum phase transition from 
a gapped Z2 spin liquid state to a gapless state with long- 
range magnetic order. 1415 The slave particles can also be 
fermions 16-19 . Hastings 16 classified possible states with bro- 
ken symmetries, including valence-bond solids and liquids 
starting from a parent state with relativistic spinons as low 
energy excitations. Later works 17,18,20 concluded that the par- 
ent state is the ground state by optimizing numerically the 
Gutzwiller-projected mean-field solutions using variational 
Monte Carlo. Based on an projective symmetry group anal- 
ysis, Lu et al [9 suggest that one of the candidate ground states 
is similar to the Z2 liquid state found in DMRG calculations. 
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It is clear that the choice of particle statistics at the starting 
point (1) determines the statistics of spin excitations, at least 
at the level of the mean-field approximation. It is therefore de- 
sirable to determine the statistics of excitations at the outset by 
some other means. In our previous work, 21 we presented ar- 
guments that the Heisenberg antiferromagnet on kagome has 
fermionic spinons. To do so, we constructed states with two 
spin- 1/2 solitons that resemble closely spinons of the one- 
dimensional antiferromagnet on a sawtooth chain 22,23 . We 
showed that the wavefunction is antisymmetric under the ex- 
change of two spinons, hence fermionic statistics. This calcu- 
lation was performed on a tree version of the kagome lattice. 24 
The deformation of the lattice was necessary to control the 
number of spinons. Strong attraction, mediated by exchange 
interaction, binds two spinons into a small pair with spin zero, 
which manifests itself as a defect — a triangle lacking a va- 
lence bond. A kagome lattice proper has a finite concentration 
of such defects: one in four triangles. This counting yields one 
spinon per unit cell on kagome. The main goal of this paper 
is to build a many-body model of interacting spinons of this 
kind. 

Our approach is variational in essence. One part of the 
variational basis consists of "dimer-covering" states where all 
spins form singlets, "dimers", with one of their nearest neigh- 
bors. The wave function of such a dimer between site i and 
site j is: 

|ii>6 = ^(|it,i;>-|it,H». (2) 

To prevent a sign ambiguity, the dimer state \ij)b is repre- 
sented with an arrow pointing from site i to site j. This is the 
bosonic convention 25,26 for dimer wave function as we can 
write it by using the Schwinger bosons to represent spins: 

\ij) b = 2-^(bl t bl-b\/ n )\0) = -\ji) b (3) 

The reason to includes such states is as follows. The Heisen- 
berg Hamiltonian on a kagome lattice can be rewritten in 
terms of the total spins of individual triangles Sa- 

H = J^S i .S j = ^j:(sl-l). (4) 

(ij) A V 4 / 

For spins of length 5=1/2 the energy would be minimized 
if every triangle had the lowest possible spin, 5a = 1/2. This 
can be achieved if two of the three spins on every triangle 
form a singlet (a quantum dimer), leaving the third spin to 
form a singlet on another triangle, Fig. 1. This program can 
be realized for a one-dimensional analog of kagome, the saw- 
tooth chain, 22 ' 23 which has two dimerized ground states. Un- 
fortunately, the trick does not work on kagome, where one 
quarter of triangles lack a dimer. 1 In our previous work, 21 
we demonstrated that these defect triangles are bound states 
of two fermionic spinons with total spin 0. This translates 
into two spinons for every four triangles, or one spinon per 
unit cell of kagome. Although this is different from the stan- 
dard slave-fermion approach, in which there are three spinons 
per unit cell (one per site), the two pictures are closely re- 
lated as we discuss in later sections. Our spinons live on the 



honeycomb lattice formed by the centers of kagome triangles. 
Their motion is strongly constrained by the presence of quan- 
tum dimers. These constraints can be described in terms of 
an emergent compact U(l) gauge field. In addition to inter- 
acting with the gauge field in the usual manner, our spinons 
experience strong short-range attraction to each other in the 
5 = channel. This is captured by an on-site negative Hub- 
bard U interaction between spinons. Because attraction be- 
tween fermions generally induces Cooper pairing, the low en- 
ergy theory of our model is a compact U(l) gauge theory in- 
teracting with a charge-2 Higgs field, which has a Z2 liquid 
ground state. 27,28 Excitations are fluxes of the Z2 gauge field 
(visons) and deconfined fermionic quasiparticles (spinons). 

The rest of the paper is organized as follows. In Sec II A, 
we briefly review the solution of two spinons on the Husimi 
Cactus, a tree of corner sharing triangles. We then review the 
arrow representation of dimer coverings and extend it to in- 
clude spinons. The model is written down in Sec II D . We 
spent the next few sections to solve the model in various lim- 
its. We first study the model in U = limit in Sec III and 
obtain two saddle point solutions: the zero and it flux phases. 
For U large enough, both saddle points host a Z2 liquid phase 
as the ground state. Since it is more natural to realize small 
spin gap in the it flux phase, we identify it as the saddle point 
that better describes the kagome. We calculate the lower edge 
of two-spinon continuum of the Z2 liquid originated from the 
7r flux phase. Based on previous works, we discuss the na- 
ture of the finite temperature phase transition to the Z2 liquid 
phase. We conclude the paper with discussion of the relation 
between our works and previous works, possible experimental 
and numerical comparisons as well as future directions. 

n. THE CONSTRUCTION OF THE MODEL 

In this section, we build the model. We first motivate the 
elements of the model by reviewing previous works. We then 
write down the model and explain its basic properties. 

A. An isolated defect triangle on the Husimi cactus 

Before we describe the many body model, it is instructive 
to briefly review the problem of an isolated defect triangle on 
the Husimi cactus 21,29 . 

The Husimi cactus (Fig 2) is the Cayley tree of corner shar- 
ing triangles. It reproduces the local geometry of kagome 
without the presence of loops. The ground state of the Hamil- 
tonian (4) is any dimer covering state since all triangles are 
vacuum triangles of lowest possible energy. An defect trian- 
gle can be isolated at the center of the Cactus. We label the 
state |0, 0, 0). 

By applying the exchange Hamiltonian, we map out the full 
Hilbert space, T2 ,that |0, 0, 0) belongs to. Generally, a state 
in T2 is characterized by the position of two mobile spin 1/2 
particles, spinons, connected by a long range singlet. The two 
spinons move on three one dimensional trails x, y and z 29 con- 
nected to the center of the cactus. Such a state is labeled as 
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FIG. 2: (a) An isolated defect triangle at the center of the Husimi 
Cactus, (b) The defect triangle is broken up into two spinons con- 
nected by a long range singlet (the blue dashed line). (c),(d) The 
exchange of two spinons reverse the "direction" of one singlet. 



\x, y, z) with the constraint that xyz = 0. The dimers con- 
strain the motion of spinons by determining the three one di- 
mensional paths. 

We project the exchange Hamiltonian (4) to the IY For 
example, we have: 

H\x, y,0)=-t(\x+l,y,0) + \x- 1, y, 0) + \x, y + 1, 0) 
+|s,y-l,0))-[tf x ,o* v ,o|0 > ) 0) 

(5) 

where t — J/2 and U = 3 J/4. The spinons tend to delocalize 
to gain kinetic energy. On the other hand, if they form a near- 
est neighbor singlet adjacent to a defect triangle, the potential 
energy of the system is lowered by U. We obtained the spec- 
trum of the two spinons with total S = 0. 21 ' 29 In addition to 
the two particle continuum, the two spinons can form a bound 
state gaining energy 0.06 J. 

The cactus also provides us with an opportunity to deter- 
mine the statistics of spinons by performing their adiabatic 
exchange. 21,29 After an exchange, the wave function acquires 
a negative sign because an odd number of dimers change from 
\ij)b to \ji)b = —\ij)b- The fermionic statistics of spinons 
comes about in a rather tortuous way in the bosonic conven- 
tion for spin singlets. It follows much more naturally if we 
use a fermionic convention 25 26 instead: 



\ij) f = 2-^(al t al 



a W Jt 



)|0) = \ji) fl 



(6) 



where creates a fermion with spin a on site i. Under this 
convention, the phase ambiguity is absent as \ij) j = \ji) f. 
Spinon motion can be described in either bosonic or fermionic 
representation of spin singlets (Appendix A). 

On kagome lattice, there is a finite concentration of spinons. 
These fermions have two opposite tendencies: delocalizing to 
lower kinetic energy and binding into nearest neighbor dimers 



to gain potential energy. Their motions are also constrained 
by the dimer configurations. To characterize the constrains, 
we introduce the arrow representation. 



B. Arrow representation 

Zeng and Elser 30 , and later Misguich et a/., 1,32 used an ar- 
row representation for dimer coverings on kagome. The ar- 
rows live on links of a dual honeycomb lattice, whose sites are 
centers of kagome triangles. When a quantum dimer is present 
on a triangle, two arrows point into this triangle through the 
ends of the dimer. Thus a triangle containing a dimer has two 
arrows pointing in and one out. A defect triangle has three ar- 
rows all pointing out. A spinon (of the antikink type 21 - 29 ) has 
one arrow in and two out, Fig. 1. 

In the arrow representation, spinons live on honeycomb 
sites and move in the direction of arrows only. Thus the spinon 
in Fig. 1 (b) can only move up or to the left, but not to the right. 
As a spinon moves across a link, the arrow on that link re- 
verses its direction. These rules preserve the right arrow count 
(one in, two out) for the moving spinon. In what follows, we 
use boldface indices r to represent sites of the honeycomb lat- 
tice and italic indices i for kagome sites. 



C. Compact U(l) gauge theory 

The arrow representation can be reparametrized as a U(l) 
gauge theory on the honeycomb lattice. We define a gauge po- 
tential A rr > = — A r > r on a link connecting neighboring hon- 
eycomb sites r and r'. The gauge field is compact in such a 
way that its wavefunction is antiperiodic: ip(A rr i + 2ir) = 
—ip(A Tr '). The momentum conjugate to the gauge field, 
the link electric field E rr > = —id/dA rT i = —E r / T , takes 
on half-integer values ±1/2, ±3/2, ±5/2,... We identify ar- 
rows on honeycomb links with the smallest (in the absolute 
sense) values of the electric field, E rT i = ±1/2. The elec- 
tric charge is defined as the electric flux emerging from a 
site, Q r = Ylr'^rr'- By this definition, a triangle with a 
quantum dimer (and thus in a ground state) carries charge 
Qo = —1/2, which can be regarded as a background charge 
of the vacuum. A honeycomb site with a spinon has charge 
Q = +1/2. A defect triangle (a spinon pair) has Q = +3/2. 
Subtracting the background charge, we find that spinons carry 
charge Q — Qo = +1 and defect triangles (spinon pairs) 
Q-Qo = +2. 

The constrained dynamics of spinons, including their at- 
traction in the 5 = channel, is captured by the following 
Hamiltonian: 




U J2 



(7) 

The first term in the Hamiltonian is a compact U(l) gauge 
theory with fermions of charge 1. The E^ T ,/2e term re- 
moves states with electric fields exceeding the minimal ab- 
solute value of 1/2. The gauge factor e lAii> increments the 
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electric field E rr i by 1 as a fermion moves from r' to r. Note 
the absence of the standard magnetic flux term cos <&, where 
the flux $ through a hexagonal plaquette is the lattice analog 
of a line integral of the gauge field § A ■ dr. 

Even in the absence of fermion self-interaction (the Hub- 
bard U term) the theory (7) is not exactly solvable. In the 
unphysical limit of large e, the Hamiltonian becomes diago- 
nal in the gauge fields A rr i , which have no dynamics of their 
own and merely set a magnetic flux for the fermions. In this 
limit the standard slave-fermion mean-field approximation be- 
comes exact and the ground state is found by finding a mag- 
netic flux background yielding the lowest kinetic energy for 
fermions. A state with well-defined fluxes, and thus gauge 
variables A rr ', is clearly unphysical: it has infinitely large 
fluctuations of the electric fields E rr > . In the physical small-e 
limit, E rr i is limited to ±1/2 and gauge variables A rr i un- 
dergo strong fluctuations. We need to write the theory in a 
more economical form in order to make further progress. 



D. Pseudospin representation 

An alternative approach is to represent the arrows on links 
as pseudospins of length s = 1/2. Using a bipartite nature of 
the honeycomb lattice, we define s* r , = s*, r = +1/2 if the 
arrow points from sublattice A to sublattice B. The number of 
spinons on site r is related to the pseudospins on the adjoining 
links: 



(8) 



with the symbolic notation (— l) r = +1 if r is on sublattice A 
and —1 if r is on sublattice B. 

In the pseudospinon representation, the Hamiltonian is 



H 



rcr^rr' r " 



H.c.) 



(9) 

The coupling constants t and U are both of the order of J, the 
only energy scale in the Heisenberg model. A crude estimate 
based on our study of spinons on a Husimi cactus yields t = 
J/2 and U = 3 J/4. The number of spinons is 1 per unit cell. 

Local constraints (8) give rise to a compact U(l) gauge 
symmetry. To make it manifest, we write down the La- 
grangian for pseudospin variables expressed in terms of angu- 
lar variables, s z = s cos 0, s ± = ssin6e ±l ^. A path integral 
for this theory can be written as 

J DA D(j> DO Da) Da cxp (i J Ldt^J, (10) 



where 



L = s (cos 6 Tr > - l)</> rr ' + i E o) r(J a ra (11a) 

(rr') r,cr 



r 

-H. 



(lib) 

(lie) 



The first line (11a) contains the standard kinetic terms for 
spins and nonrelativistic fermions; the next term (lib) en- 
forces the local constraints. The Lagrangian is invariant under 
U(l) gauge transformations, 



— > a rtJ exp(zA r ), (12a) 
-> <t> Tr > + (-l) r (A r - A r 0, (12b) 
-> A° r + X r . (12c) 



It is evident that quantities (— l) r (/> rr ' transform like spatial 
components of a lattice U(l) gauge field and A® as its time 
component. The angular nature of <j> TT i makes the gauge field 
compact. (9) and ( 1 1) are our main results. The theory will be 
solved in various limits in the remainder of the paper. 



III. THE (7 = LIMIT: GAUGE MEAN-FIELD THEORY 

We first solve the theory in the [7 = limit. It is still non- 
trivial since the fermions are interacting with a compact U(l) 
gauge field. To make progress, we decompose the kinetic en- 
ergy terms of fermions using the gauge mean-field theory: 33 



(13) 



The averages (s^ r ,) = s s'm 8 rr r exp(i<f) rr >) play the role of 
hopping matrix elements for spinons, whereas bond averages 
(a\ a a r i<j) create a transverse magnetic field for pseudospins. 
The spinon kinetic energy is minimized by setting 6 rr / = tt/2 
for every bond. The spinon spectrum is then dependent on 
azimuthal angles rr ' that determine magnetic fluxes through 
hexagonal plaquettes. Below we consider flux configurations 
that do not break the time reversal symmetry (fluxes and it) 
and preserve translational symmetry (all fluxes are the same). 



A. Zero flux phase 

In the simplest case, magnetic fluxes on all plaquettes are 
zero. The mean-field Hamiltonian for the fermions is: 

Hmf = -ts ^2 ( a ta a r'cr + H. c.) - fi^2 al a a r<T . (14) 

(rr'), a r,a 

The chemical potential fj, is adjusted so that the number of 
fermions is 1 per unit cell. The primitive vectors are ai = x 
and &2 = x/2 + s/Zyjl. After a Fourier transformation, 



H = £ Cr(k)cj C(TT c k[Tr 



(15) 



k,<7.T 



where 



k x y3ky 



£1,2 (k) = TtS\ 3 + 2 cos k x + 4 cos — cos 



(16) 
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FIG. 3: The spinon Fermi surface (black) for zero-flux phase. The 
red-dashed line is obtained by assuming E = fcj /2m. 



FIG. 5: The positions of the Fermi points pi,2- The black hexagon 
is the first Brillouin zone of the honeycomb lattice. The purple rect- 
angle is the Brillouin zone of the 7r-flux state. Triangles mark Dirac 
points pi and p2. 
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FIG. 4: M rT r — —1 for the dashed bonds. 



The chemical potential [i sw — 0.74i yields 1 fermion per unit 
cell. Spinons have a nearly circular Fermi surface (Fig. 3) 
with the Fermi momentum kf ~ 2.7. The effective mass is 
m w 4.7/t and the density of states is m/2ir. 

Thanks to the extended Fermi surface, the spinons remain 
deconfined 34 even if the compact U(l) gauge fluctuations are 
taken into account. We do note that a recent report 35 claims 
the contrary that compact gauge fluctuations confine spinons 
with a Fermi surface. 



Here M rr < = ±1 for solid (dashed) bonds in Fig. 4. The 
resulting bands are, in the order of increasing energy, 

6,2 (k) = -tsy3 ± 2/>(k) - (18a) 
&,4(k) = +tS^3 T 2p(k) - fx, (18b) 



where p(k) = y 1 + cos 2 fc x — sin fc^ sin \/3k y . The chemi- 
cal potential /.t = — 73i/2 is set at the intersection of bands 
£i and £2, so that only the lowest band £1 is filled. The 
lower two bands touch at two Dirac points with momenta 
P1.2 = ± (7r/2, 71-/273) (Fig 6). In the vicinity of the two 
Dirac points P1.2, the Fermi velocity is t/2^/2. 

The "flavor" number of Dirac fermions is Nf = 4 con- 
sidering both the spin and valley degeneracy. In the limit of 
large Nf i6 , Dirac fermions are not confined by compact U(l) 
gauge field. There are conflict reports 37 ^* on whether Nf = 4 
is above or below the critical number flavors. We will oper- 
ate under the assumption that spinons are deconfined. Since 
spinons eventually form Cooper pairs, this assumption is not 
essential to our conclusion. 



B. The 7r flux phase 

The unit cell is doubled to include four sites of the honey- 
comb lattice. The mean-field Hamiltonian for spinons reads: 

Hmf = ~ts ^2 M rr ,(al a a r , a + H. c.) - a^cw 

(rr').a r.er 

(17) 



IV. [I/O: BCS AND THE HIGGS PHASE 

We consider now the case of finite U. For sufficiently 
strong attraction, fermions form Cooper pairs and develop a 
gap around the Fermi surface. We calculate the BCS gap for 
both 0- and 7r-flux states. 
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FIG. 6: Energy bands of noninteracting spinons with the wavenum- 
ber k = (kcos<f>, fcsin^) with = n/6 in the 7r-fiux phase. Left 
panel: our model, Eq. (18). Right panel: slave-fermion theory, 
Eq. (B4); the flat band is doubly degenerate. Filled (empty) bands 
are shown by solid (dashed) lines. 



A. Zero-flux phase 

The interaction term is decoupled in the usual BCS way: 

2NA 2 

--A 

k.r 



U 



■-A$>U c -kr|+H.c.), (19) 



where N is the number of unit cells and r = 1, 2 is the band 
index. 

A Bogoliubov transformation can be used to diagonalize 
the Hamiltonian: 



H = ^e r (k) c^jAt, 



(20) 



where e T (k) = y^|(k) + A 2 is the excitation energy of 
a Bogoliubov quasiparticle. The self-consistency condition, 
A = U(a r ^a r i), yields the gap equation: 



-E 

IN t-~t 



u 



(21) 



The presence of extended Fermi surface implies that A ^ 
for infinitesimal U . For the bare couplings, t = J/2 and 
U = 3 J/4, we obtain A = 0.39t. The spin gap is then 2A = 
0.39J, much larger than the DMRG value 0.1J. In order to 
reduce the spin-gap to numerically observed value, large mod- 
ifications of U or t from their bare values are needed, which 
is not likely. 



B. The 7r-flux phase 

Similar calculations can be carried out for the 7r-flux phase. 
The gap equation reads: 



-E 

AN ^ 



V 



4N Itr v/£r(k) 2 +A2 



(22) 



Due to vanishing density of state around the Dirac points, the 
gap equation is only satisfied for U > U c where 



i = -E 

AN 



47V^|£ T (k)| 



(23) 
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FIG. 7: The lower edge of the two particle continuum for the 7r-flux 
phase assuming A = O.lt. 



Explicit calculations show U c « t. For the bare couplings, 
t = J/2 and U = 3J/4, the spin gap 2A w 0.32J. While 
it is still large, relatively small adjustments to parameters can 
improve the agreement greatly. For example, for an enhanced 
spinon hopping t = 0.68 J and the bare U = 3J/4, we find 
the spin gap 2 A = 0.1 J. 

Figure 7 shows the lower edge of the two-spinon continuum 
for A = O.lt in the 7r-flux phase along high-symmetry direc- 
tions. Our theory predicts low energy S = 1 excitations at 
both T and M points (Fig 5). This can be understood without 
explicit calculation. The low energy spinon excitations are of 
momentum pi,2- Combining with suitable reciprocal lattice 
vectors, the total momentum of two spinons with momentum 
Pi and pj correspond to momentum of V and M points. It 
should be borne in mind that the scattering intensity will be 
reduced by the antiferromagnetic structure factor (recall that 
our spinons are coarse-grained degrees of freedom). 

To understand the nature of the "superconducting" phase, 
we integrate out the fermions. If we neglect the amplitude 
fluctuations of A, the low energy theory can be described as 
compact U(\) gauge field interacting with a charge 2 Higgs 
field. The Euclidean action for the theory in d = 2 + 1 is 27,28 

S = /3^cos#(r') + hJ2 C0S ( S ^( r ) - QM r ))- (24) 

r' r.fi 

Here the summation of r' goes over all plaquettes. 8(r') is the 
flux through plaquette r'. <fi(r) is the phase of the Higgs field 
at site r. Q = 2 is the charge of Higgs field. 6* M (r) is </> r ,r+/j 
in our language. 5^ is the difference operator defined as: 



^/(r) = (-l) r (/(r + r r A)-/(r)). 



(25) 



A Monte-Carlo study 28 of model (24) determined its phase 
diagram. At low temperature (j3, h — > oo with j3/h fixed) the 
theory enters a Higgs phase which preserves local Z 2 gauge 
symmetry. In other words, the ground state of the theory is a 
Z2 liquid. 
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We estimate (3/h ~ t 2 /A 2 1. Based on the phase dia- 
gram (Figure 5 of Bhanot et a/. 28 ), we expect the finite tem- 
perature phase transition into the Z 2 liquid phase belongs to 
the three-dimensional XY universality class. The transition 
temperature is determined by the spin gap. The conclusion 
does not depend on our saddle-point choice of 7r-fiux phase. 
It is only based on that t ~ J is much larger than the spin gap. 



V. U > t LIMIT: LIMITATIONS OF THE MODEL 

In the U —> 00 limit, all dimer covering states have the 
same energy. Using these states as our basis, we employ the 
degenerate perturbation theory to calculate the energy contri- 
bution of the hopping term to the lowest order in t, i 6 /[/ 5 in 
this case. 

Without explicitly calculating these terms, their general fea- 
tures already reveal the limitation of our model. These opera- 
tors shift dimers around a star collectively. The amplitudes of 
the shifts is: 



compounds, the presence of the DM interaction and of mag- 
netic impurities will affect the two-spinon edge. While com- 
paring with future measurements on more "ideal" materials is 
certainly desired, it is crucial to cany out calculations taking 
into account such perturbations, which will be the focus of the 
future work. 
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Appendix A: The motion of a spinon on the sawtooth chain 

We consider a spinon on an one-dimensional chain of cor- 
ner sharing triangles, the saw tooth chain 22,23 . Its equation of 
motion is derived using both bosonic and fermionic conven- 
tions of dimer wave functions. 



(26) 



1. The bosonic convention 



where c(h) > is the symmetry factor depending on the num- 
ber of dimers and their configurations. 

These amplitudes have been derived from S = 1/2 Heisen- 
berg antiferromagnetic model using the overlap expansion. 
Our results do not match such calculations. In other words, 
our model is not expected to give a good descriptions of sin- 
glet excitations. 



VI. DISCUSSION 

We first comment on connections between our model and 
the slave-fermion theories, 16-19 particularly the algebraic spin 
liquid. 1718 Like the previous authors, we focus on a phase 
with a U(l ) flux ir on hexagonal plaquettes and find, for a sim- 
ilar gauge choice, low-energy fermion excitations near mo- 
menta P12 = ±(7r/2,7i"/2\/3). Furthermore, the four bands 
in our model have the same dispersions (aside from an overall 
scale factor) as the four upper bands 16 of the algebraic spin 
liquid (Appendix B). In this sense, it seems that our theory 
is an economical low-energy description of the slave fermion 
theory. The additional two bands in the slave-fermion models 
are dispersionless and could be related to spinons of another 
flavor, 21 which also happen to be dispersionless in the pure 
Heisenberg model. 

In our approach, the compact U(l) gauge theory emerges as 
a natural way to incorporate the constraint of dimers on spinon 
motion. In contrast, the gauge field is used to enforce the 
constraint on spin length, which translates into a requirement 
of one spinon per site. All states in our variational basis satisfy 
this constraint. 

Our model shows that low-energy triplet excitations exist at 
both T and M points. It also predicts the universality class of 
the quantum phase transition into the Z 2 liquid phase. In real 



We define the state with a spinon on site 2n to be: 

\n,*) b = &L.J0) n i 2 i' 2 i+ i > n i 2 ^ 2 -? - !) < ai > 



j>n 



where b\ n a creates a bosonic spinon of spin a at site 2n. 
The exchange interaction for bond (ij) can be written as: 



13 2 3 4 



(A2) 



where Pjj permutes the spin state of site i and j. Applying 

H 2n ,2n+i to \n, a) b , we get: 

H 2n ,27i+i\n) b = --^\ n + l)b + j\ n )b- (A3) 

The full equation of motion reads 21 ' 29 : 

7 5 / 

H\n) b = --(\n + l) b + \n- 1)„) + —\n) b . (A4) 



2. The fermonic convention 

We again define: 

l». *)f - 4n,J0) J] \2j, 2j + l) r J] \2j, 2j - l) f . (A5) 

j<n j>n 

The full equation of motion can be derived similarly: 

7 5 7 

H\n)j = g (l»» + !)/ + \ n ~ %) + ~\ n )f- ( A6 > 

The spinon hopping amplitude changes sign comparing to 
equation A4. Similar sign-change persists on both the Husimi 
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Cactus and the kagome lattice. In our model, spinons live on 
the honeycomb lattice with plaquettes of even length. The 
overall sign of spinon hopping amplitude is thus irrelevant 
since it can be reversed by a simple gauge transformation. 

Appendix B: Slave-fermion mean-field theory on kagome 

In this section, we solve for the mean-field band structure 
of the algebraic spin liquid. In the slave fermion theory, a spin 
operator S^ a ' is written in terms of fermions: 

s { t a) = (Bl) 

where with a = x, y, z are the Pauli matrices. 
The exchange Hamiltonian can be written as: 

H = - 4a a i"' a ]a' a jv ( B2 > 

We decompose the four-fermion terms using mean-field 



approximation : 



hmf ~ ^Bwi^ +ILc -) + 2 5 i^i 2 ' (B3) 

(ij) (ij) 



with the self-consistency condition Xij = ~( a ia a j<r)- By 
symmetry, \xij\ = X f° r every nearest-neighbor bond (ij). 
Assuming that fluxes through hexagons and triangles are it 
and 0, respectively, we obtain six bands, in the order of in- 
creasing energy: 

a, 2 (k) = -J X , (B4a) 
6,4k) = ^(1-V3±2^), (B4b) 

&, 6 (k) = ^(1 + v/3T2p k ). (B4c) 
These are shown in the right panel of Fig. 6. 
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